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Abstract 

We prove, for the relativistic Boltzmann equation on a Bianchi Type 
I space-time, a global existence and uniqueness theorem, for arbitrarily 
large initial data. 


1 Introduction 

The Boltzmann equation is one of the basic equations of the kinetic theory. This 
equation rules the dynamics of a kind of particles subject to mutual collisions, 
by determining their distribution function, which is a non-negative real-valued 
function of both the position and the momentum of the particles, the distribu¬ 
tion function, also called the ’’phase space density” of the particles, is physically 
interpreted as ’’the probability of the presence density” of the particles in a 
given volume, during their collisinal evolution, and is an essential tool for such 
a statistical description. 

In the case of instantaneous, localized, binary and elastic collisions, the 
distribution function is determined by the Boltzmann equation, through a non¬ 
linear operator called the ’’collision operator”, that acts only on the momentum 
of the particles, and that describes, at any time, at each point where two parti¬ 
cles collide with each other, the effects of the behavior imposed by the collision 
to the distribution function, taking into account the fact that the momentum 
of each particle is not the same, before and after the collision, only the sum of 
their two momenta being preserved. 

Due to its importance in the kinetic theory, several authors studied and 
proved local and global in time existence theorems for the Boltzmann equation, 
in both the non-relativistic case, that considers particles with low velocities, 
and the full-relativistic case, which includes the case of fast moving particles 
with arbitrarily high velocities, such as, for example, particles of ionized gas 
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in some medias at very high temperature as : burning reactors, solar winds, 
nebular galaxies. 

1) In the non-relativistic case, the original global result is due to Carleman, 
T. , in UJ ; Diperna, R. J. and Lions, P. L. proved global existence and 
weak stability in (3j ■ Di Blasio, G. , proved the differentiability of spatially 
Homogeneous Solutions of the Boltzmann Equation in the non Maxwellian 
case in Illner, R. , and Shinbrot, M. , proved a global result in 

in the case of small initial data and without symmetry assumption ; an 
analogous result is unknown in the full relativistic case. 

2) In the full relativistic case, several authors proved local existence theorems, 
considering this equation alone, as Bichteler, K. , in [5], Bancel, D. , in 
PQ, or coupling it to other fields equations as Bancel, D. , and Choquet- 
Bruhat, Y. , in j2]. Glassey, R. , T. , and Strauss, W. , obtained a global 
result in (Jj , in the case of data near to that of an equilibrium solution with 
non-zero density. Noutchegueme, N. and Tetsadjio, E. M, proved Global 
existence for small initial data on the Minkowski space-time in m 

The objective of this paper is to extend to the full relativistic case, the 
global existence theorem in the case of arbitraliry large initial data. That 
was certainly one of the goals of Mucha, P. , B. , who studies in m and 
GU, the relativistic Boltzmann equation, coupled to Einstein’s equations, 
for a flat Robertson-Walker space-time. Unfortunately, several points in 
that work are far from clear, such as : using a formulation which is valid 
only in the non-relativistic case . 

In this paper, we begin, first of all, by giving, following Glassey, R. , T. , in 0, 
the correct formulation of the relativistic Boltzmann equation, taking as back¬ 
ground space-time a Bianclii type I space-time . We consider the homogeneous 
case, which means that the distribution function only depends on the time 
and the momentum of the particles. Such cases are very useful for instance in 
cosmology. Next, we build a suitable functional frame-work, and we construct 
a sequence of operators to approximate the ’’collision operator”; this provides 
a sequence of real-valued functions that converges to a global solution of the 
Boltzmann equation. We prove this way, using very simple a priori estimates 
that lead to a considerable simplification of the method followed in M, that, 
if the coefficients of the back-ground metric are bounded away from zero, then 
the initial value problem for the relativistic Boltzmann equation on that Bianchi 
type I space-time, has a unique global solution, for arbitraliry large initial data 
and that the solution admits a very simple estimation by the initial data. 

The paper is organized as follows : In section 2, we specify the notations, 
we define the function spaces, we introduce the relativistic Boltzmann equation 
and the ’’collision operator”. We end this section by a sketch of the strategy 
adopted. In section 3, we give the properties of the ’’collision operator”, some 
preliminary results, and we construct the approximating operators. Section 4 is 
devoted to the global existence theorem. 
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2 Notations, function spaces and the relativistic 
Boltzmann equation 

2.1 Notations and function spaces 

A Greek index varies from 0 to 3 and a latin index from 1 to 3 . We consider 
as back-ground space-time, a Bianchi type I space-time denoted (R 4 , p), where, 
for x = ( x a ) = (x°, x ®) £ R 4 , x° = t, denotes the time and x = (x l ) the space ; 
g stands for the metric tensor with signature (—, +, +, +) that can be written : 

g = — dt 2 + a 2 (t)(di 1 ) 2 + b 2 (t)[{dx 2 ) 2 + (da; 3 ) 2 ] (2.1) 

in which, a(t) and b(t) are given non-negative regular, real-valued functions. We 
require that a(t) and b(t) are bounded from below i.e. : 

a(t) > C 0 , b(t) > C 0 (2.2) 

where C a > 0 is a given constant. Let us obseve that for a = b, (R 4 ,g) reduces 
to the fat Robertson-Walker space-time, which is the basic model for the study 
of the expanding Universe. We consider the collisional evolution of a kind of 
uncharged particles in the time-oriented curved space-time (R 4 ,g). An essential 
tool to describe the dynamics of such particles is their distribution function we 
denote by /, and that is a non-negative real-valued function of both the position 
(x a ), the 4-monrentum p = ( p a ) = (p°,p z ) of the particles, and that coordinalize 
the tangent bundle T(R 4 ) i.e : 

/ : T(R 4 ) ~ R 4 x R 4 -> R+, ( x a ,p a ) i—* f(x a ,p a ) £ R+. 

We define a scalar (or dot) product and a norm on R 3 , by setting, for 
p = (p°,pO = (p°,p) and q = ( q °, q l ) = {q °, q) : 

p.q = a 2 (t)p Y + b 2 (t)(p 2 q 2 +p 3 q 3 y, \ p | 2 = a 2 (t)(p 4 ) 2 + 6 2 (t)[(p 2 ) 2 + (p 3 ) 2 ] 

(2.3) 

We consider massive particles with a rest mass that can be rescaled to m = 1. 
The particles are then required to move on the future sheet of the mass-shell 
whose equation is g(p,p) = — 1, or equivalently, using (12.111 and 12.31) : 

p° = Vl+|p| 2 (2.4) 

The relation itlOl) shows that / is in fact defined on the subbundle of T(R 4 ) 
whose local coordinates are x a and p L . 

Now we consider the homogeneous case for which / depends only on the time 
x° = t and p = {p l ). The framework we will refer to is L 4 (R 3 ) whose norm is 
denoted || . || ; we set, for r £ R, r > 0 : 

Xr = {f£L\ R 3 ), /> 0, me., 11/|| < r}. 

Endowed with the distance induced by the norm || . ||, X r is a complete and 
connected metric space. Let I be a real-interval; we denote by C[I; L 1 (R 3 )] the 
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Banach space of continuous and bounded functions g : 1 i—> L 1 (R 3 ), endowed 

with the norm ||| g |||= Sup || g{t) ||. We set : 
tel 

C[I ; X r ] = { g G C[I; ^(R 3 )], g(t) £ X r , Vi € /} 

Endowed with the distance induced by the norm ||| . |||, C[I;X r ] is a complete 
metric space. 


2.2 The relativistic Boltzmann equation 

The relativistic Boltzmann equation for uncharged particles, in the homogeneous 
case, on the Bianchi type I space-time (R 4 ,g) we consider , can be written : 


dt a dp 1 





3 df 1 




(2.5) 


with a = 4|, and where Q is the collision operator we now introduce. In the 
case of instantaneous, localized, binary and elastic collisions we consider, the 
collision operator Q that acts only on the momentum variable, is defined as 
follows, p and q standing for the momenta of 2 particles before their collision, 
p' and q' standing for their momenta after the collision, regardless for the time 
i, and where / and g are 2 functions on R 3 : 


i) Q(f,g)=Q + (f,g)-Q-(f,g) (2.6) 

where 

2) Q + (Lg)=[ ^ dq f f{p')g{q')S{a(t),b{-t),V,cLp'iq')du 

J R3 q° Js 2 

(2.7) 


3) Q (f,g) = ( a ^ b ^ dq I f(p)g(q)S(a(t), b(t),p, g,p', q')dio 

J rs q Js 2 

( 2 . 8 ) 

in which 

4) S 2 is the unit sphere of R 3 whose area element is denoted dcu ; 

5) S is a non-negative real-valued function of the indicated arguments, called 
the collision kernel or the cross section of the collisions, and for which we require 
the boundedness and the symmetry assumptions : 

0< S(a(t), b(t),p, q,p f , q') < Ci (2.9) 


S(a(t), b(t),p, q,p', q') = S(a(t,),b{t),p\q\p,q) ( 2 . 10 ) 

where C\ > 0 is a given constant. 

6) As consequences of the conservation law p + q = p' + q\ that gives : 

P° + q° = (PT + (q')° ; p + q=p' + q‘' (2.11) 
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i) We have, using (2.4) and the first equation (2.11) : 


+ Vi+TW = V^+Wf + V^+Wf (2.i2) 

which expresses the conservation of the quantity : 

e = y/l+ I p I 2 + \/l+ I q I 2 (2.13) 


called the energy of the unit rest mass particles. 


ii) We set, to express the second equation and following Glassey, R. 

T., in |8J : 

Jp' =P+C(p, q,u)w 
1 q' = q-C(p, q,u)u> 


(2.14) 


in which C(p , q , u>) is a real-valued function . The relations 12.1211 and 
mu lead to a quadratic equation in C that solves to give the only non¬ 
trivial solution : 


C(p,q,w) = 


2 p°q°e uj.(q — p) 
e 2 - [u>.(p+ q)] 2 


(2.15) 


with p = -2;, e defined by 12.1311 and the dot product by (IQ) . Now, 
using the usual properties of determinants, the Jacobian of the change of 
variables (p, q) i—> (p', q') defined by 12.1411 is computed to be : 


d{p',q') _ p'°q'° 

d{p , q) P°q° 


(2.16) 


Notice that formulae 12.1511 and trna are generalizations to the considered 
Bianchi type I space-time, of analogous formulae established in Glassey, R. 
T., {5J, in the case of the Minkowski space-time, whose metric corresponds 
to the particular case of mu, when a[t) = b(t) = 1. 


It then appears clearly, using mi and mu , that the functions in the inte¬ 
grals ma and mi depend only on p, q , u> and that these integrals give func¬ 
tions Q + (f,g), Q~{f,g), of the single variable p. In practice, we will consider 
functions / on R 4 that induce for t £ R, functions f(t) on R 3 , defined by : 

f(t)(p) = f{t,p). 

Now solving the Boltzmann equation urn which is a first order p.d.e. , is 
equivalent to solving the associated characteristic system, that writes, taking t 
as parameter : 


dp 1 

dt 


2 a dp>_ 

a dt 



dj_ 

dt 


-Q(fJ) 

p° 


(2.17) 


We study the initial value problem for the system mu , with the initial data : 


p i {0) = y i -, m = fo 


(2.18) 
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In 12.171) . the equations in p 1 integrate at once to give, setting y = (j/®) € R 3 : 

p(t,y) = A(t)y, with: A(t) = Diag ( 2 - 19 ) 

whereas the initial value problem for / is equivalent to the integral equation : 

f(t,y) = fo(y)+ [ 4 Q(fJ)(s,y)ds ( 2 . 20 ) 

Jo P 

Finally, solving the initial value problem for the relativistic Boltzmann equation 
12.511 . with the initial data f 0 , is equivalent to solving the integral equation font 
in /. The strategy we adopt to solve will consist of : 

a) Constructing approximating operators Q n , with suitable properties, that 
will converge in the L 1 (R 3 )-norm to the operator ^ Q. 

b) Solving by usual methods, the following approximating integral equation: 


f(t,y) = fo(y) + f Qn{f,f){s,y)ds (2.21) 

JO 

obtained by replacing in CEB, Q by Q n , and to obtain a global solution 
f n that will converge to a global solution / of csza , in a suitable function 
space. Notice that, in R 3 , by J2T9J, the volume elements dp, dy, and the 
corresponding norms || . ||p, and || . \\ y we adopt to be || . || , are linked by 


dp 


a 2 (0)6 4 (0) 

a 2 (f)6 4 (f) 


a 2 (0)6 4 (0) 

a 2 (t)6 4 (t) 


( 2 . 22 ) 


3 Preliminary results and approximating 
operators. 


In this section, we give properties of the collision operator Q, we construct and 
we give the properties of a sequence of approximating operators Q n to — Q. 


Proposition 3.1 : 

Iff,g e i 4 (M 3 ), then Q + (f,g),Q~(f,g) e L 1 (R 3 ), and; 


pU 

1 


- 0 Q + (.Lg) II 

< 

C || /|| || 9 || 

(3.1) 

p° 

- 0 Q~{f,g) II 

< 

C || /|| || 9 || 

(3.2) 

p° 

4 Q + (9,9 ) II 

< 

C(\\f\\ + \\g 

\\)\\f-g\\ (3.3) 

p° 

4 Q~(.9,9) II 

< 

C(\\f\\ + \\g 

4? 

CO 

1 

p° 

- 4 Q(9,g) II 

p° 

< 

C(\\f\\ + \\g 

S' 

CO 

1 
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where C = 8nCia 2 (0)b 4 (Q)/C% 

Proof 

1) The expression 12.711 for Q + (f,g) gives, using the upper bound (I2.hll of S 

II ^Q + (f,9) ||p< C ia {t)b 2 (t) [ j 44 f | f{p') | | g(q r ) | duj (a) 
P° J R3 Jr 3 p°q° J s 2 

Now, the change of variables (p,q) i—> ( p',q' ), defined by (12.141) and its 
Jacobian (l2~TTil) give dpdq = — plo q qlo dp' dq'\ then we have from (a), using 
(IPl that gives ^ p ,)J^ qlo ^ < I'¬ 
ll 4 Q + (f,9 ) ||p< Cia(t)b 2 (t) f || f{p') || dp' [ || g(q') || dq' f dui = M 

P j R3 J R3 J S 2 

where 

M = 47rCia(t)6 2 (t) || / ||p || g ||p 

and (3.1) then follows from the second relation (I2.22H and the lower bound 
(I2.2l) for a and b □ 

2) The expression (2.8) for Q~(f,g) gives, using (12.hi and (12.41) that gives 

-±-z < 1: 

pO q o _ 

II 4 Q"(/.s) Up- Cia(t)6 2 (t) [ || f(p) || dp [ || g(q) || dq f dw = M 

P J R3 J R3 Js 2 

and (3.2) then follows as above. □ 

3) Inequalities (3.3), (3.4), (3.5) are direct consequences of (3.1), (3.2), 

Q = Q + — Q ~, and the fact that, the expressions and (2.8) show 
that Q + and Q~ are bilinear operators, and so, we can write, P standing 
for j- Q+ or ± Q~ : 

P(f, f ) - P{3,9) = P{f, f ~ 9)+ P(f~ 9,9) 

This completes the proof of Proposition 3.1 □ 

We now state and prove the following result on which relies the construction of 
the approximating operators; C denotes the constant defined in Proposition 3.1 


Proposition 3.2 Let r be an arbitrary strictly positive real number. Then, 
there exists an integer n 0 (r ) > 0 such that, for every integer n > n a (r) and for 
every v £ X r , the equation : 

\fnu — —Q(u,u) = v (3-6) 

has a unique solution u n £ X r 
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Proof 

Let r > 0 be given and let v £ X r . We can write (3.6), using Q = Q + — Q~ : 
\fnu — — Q + (u,u) + — Q~(u,u) = v 

71 ° 71 ° 


Notice that one deduces from the expression (2.8) of Q~ , since p is not a variable 
for that integral, that Q~(f,g) = fQ~(l,g). Hence the above equation can be 
written, using Q~(u,u ) = uQ~(l,u), for u > 0 a.e. 

u = F n {u ) (a) 


where 


F n {h) 


V + jsQ + {h,h) 
Vn + ±Q~(l,h) 


(b) 


(a) shows that a solution u of (3.6) can be considered as a Hxed point of the map 
h i—> F n (h) defined by (b). Let us shows that F n is a contracting map from X r 
to X r . We will then solve the problem by applying the fixed point theorem . 


i) We deduce from (3.1) that F n : X r i—> L 1 (K. 3 ) , and that since v, h £ X., 

" v I 


K(h) || < 


c II h II II h II < r(l + Cr) 


(c) 


ii) Let g, h £ X r . (b) gives : 


v + -^Q + (h,h) v + -^Q + (g,g) 
F n (h)-F n (g) = - p ^ J - p 


Vn + jzQ (1 ,h) y/n+^Q (1,5) 
Since g > 0, h > 0 a.e., we have , taken n > 1: 

Wn+ ^Q~{l,h)][^/n + -^Q~(l,g)] >n>y/n.. 

P° P° 


(d) 


We then deduce from (d), using once more Q (f,g) = fQ (1,5), that : 

Vn. || F n (h) — F n {g) || <|| -^Q + {h,h) - ^Q + {g,g) || + || ^ Q~(v,g ) - ^Q~(v,h) || 

+ 11 ^Q~l^Q + ( h ’ h )’9]- ^Q-[^Q + (g,g),h] || (e) 


We can write, using the bilinearity of Q : 

-^Q~{v,g)~ -4 Q~{v,h) = ^-Q~{v,g-h) 

pO rpO rpO 

So, in (e) we apply (3.3) to the first term, (3.2) to the second term to obtain, 
using || g ||< r, || h ||< r, || v ||< r : 

|| 4 Q + (h,h )— 4 Q + (di 5) II + II -^Q~{v,g) + -^Q~(v, h ) ||< 30 || h-g (f) 

pO rpO pO pO 
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Concerning the last term in (e), we can write, using the bilinearity of Q 


^Q + (h, h),g 

p° 

-- 0 Q~ 

p° 

-^Q + (g,g),h 

p° 

p° 

■^Q + (h,h),g - h 

p° 




i 

pO 

Q- 

-Q+(h,h)~ 

p° 


We then apply to the right hand side of the above equality : 

1) For the first term, property (3.2) of Q~ , followed by property (3.1) of Q + 

2) For the second term, property (3.2) of Q~, followed by property (3.3) of 
Q + to obtain, using || g ||< r, || h ||< r : 


w 

p° 


\q + {H, h),g 

p° 


p° 


-^Q + (g,g),h 

p° 


< 3C 2 r 2 || h-g 


thus (e) gives, using (f) and (h) : 


F n (h) - F n (g) ||< 3CV(1 + Cr) || h-g 


(h) 

(i) 


It is then easy , using (c) and (i), to choose an interger n 0 (r) such that, for 
n > n 0 (r) , we have: 


1 + Cr 

y/n. 


< 1 


and 


3CV(1 + Cr) 1 
Vn “2 


(3.7) 


This shows , using once more (c) and (i) that, for n > no(r), F n is a contracting 
map of the complete metric space X r into itself. By the fixed point theorem, 
F n as a unique fixed point u n G X r that is the unique solution of (3.6) □ 

In all the following, n D (r) stands for the integer defined above. 

The above result leads us to the following definition : 

Definition 3.1 Let n £ N, n> n 0 (r) be given. 

1) Define the operator 


R(n, Q) : X r —> X r , u i—> R(n , Q)u 
as follows : for u € X r , R(n , Q)u is the unique element of X r such that : 

\/nR(n, Q)u - Q [I?(n, Q)u, R(n, Q)u] = u (3-8) 

p° 

2) Define the operator Q n on X r by : 

Q n (u,u) = ny/nR{n,Q)u — nu (3-9) 


We give the properties of the operators R(n , Q) and Q n . 


Q + (g,g),h 
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Proposition 3.3 : Let r £ R+, n,m £ N m > n 0 [r), n > n 0 {r), be given. Then 
we have, for every u,v £ X r 

(3.10) 

(3.11) 

(3.12) 


(3.14) 

(3.15) 

(3.16) 


| y/nR(n, Q)u | 

= 

II«II 

\/nR(n , Q)u — u | 

< 

K 


n 

| R{n,Q)u — R(n,Q)v || 

< 

-7= II U-V || 

\tn 

Q n {u,u) 

= 

^Q [VnR{n,Q)‘ 

II Q n {u,u) - Qn(v,v) | 

< 

K || u — v || 



K 

II Qn{u,U )- -Q(v,v) || 

< 

-1 - K \\ u — u 

/pU 


n 

|| Qm(v i v') || 

< 

K | u - v || + — 
n 


Here K = iv (r) is a continuous function of r. 

Proof. 

1) (3.10) will be a consequence of the relation: 


[ ^Q(f,g)dp = 0, Vf,g£L\R 3 ) (3.17) 

J R3 P° 

we now establish. This is where we need assumption (2.10) on the kernel S. 
We have, using definition (2.7) of Q + , the change of variables (p,q) '—> (jJq 1 ) 
defined by (2.14) and its Jacobian (2.16), and since p° = y/l+ \ p | 2 : 


I n = 


f 1 

I R3 P° 


Q + (f,g)dp = a{t)b 2 (t) / dp' / dq’ 


f(p')9{q')S{p',q',p, q) 


I R3 


Is? 


s 2 v /1 + I p 1 I 2 \A+ I q' I 2 

On the other hand, we have, using definition (2.8) of Q~ and p° = \/l+ \ p \ 2 : 


du 


Jo = 


/ fQ-u, S )dp=a(tm)[ dpt dpt 

/R3 P J R3 J R3 Js 2 V 1+ I P I 2 V 1+ | q I' 




It then appears that I 0 = J 0 and 13.171 1 follows from Q = Q + — Q . Then, 
to obtain (3.10), we integrate (3.8) over IR 3 with respect to p, use (3.17), and 
conclude by 12.2211 that gives the equivalence between || . ||p and || . || □ 

2) The definition 13.811 of R(n, Q) gives, using property (3.5) of Q with 
/ = R(n , Q)u and g = 0 : 

|| V nR(n , Q)u - u || = || ^-Q [R{n, Q)u, R{n , Q)u] ||< G || R{n , Q)u || 2 

p 

and (3.11) follows from (3.10) that gives || R(n,Q)u ||= 

3) Subtracting the relations (3.8) written for u and v gives: 

u — V I 1 

R(n, Q)u—R(n , Q)v = — -= — Y—j=—Q [R(n , Q)u, R(n, Q)u] - -=—Q [R{n, Q)v , R(n, Q)v] 


yfn sfnp° 


y/np° 
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which gives, using property (3.5) of Q and once more (3.10): 



H-|| R{n, Q)u - R{n 1 Q)v || 

Jn 


(3.12) then follows from n > n 0 (r ) and (3.7) that implies ^'=■ < i. 

4) The expression (3.13) of Q n is obtained by multiplying (3.8) by n, using its 
definition (3.9) and the bilinearity of Q. 

5) Property (3.14) is obtained by using (3.13), property (3.5) of Q, followed by 
properties (3.10) and (3.12) of R{n,Q). 

6) Notice that by (3.10), sjnR(n , Q)u £ X r since u £ X r . Then, the expression 

(3.13) of Q„ and property (3.5) of Q give : 


II Qn(u , u)-^-Q(v, v) ||< 2 Cr || ^Jn.R{n , Q)u—v ||< 2 Cr (|| y/nR{n, Q)u - u || + || u — v ||) 

p° 



then, (3.16) follows from (3.5), (3.10) and (3.11), adding and subtracting u and 


v. □ 


Remark 3.1 . (3.15) shows by taking u = v, that Q n converges pointwise in 
the L 1 (M 3 )-norm to the operator XQ _ From there, the qualification of 
’’approximating operators” we give to the sequence Q n . 

We now have all the tools we need to prove the global existence theorem. 

4 The global existence theorem 

With the approximating operator Q n defined by (3.8), we will first state and 
prove a global existence theorem for the ’’approximating” integral equation: 



f(t 0 + t,y)=g 0 (y)+ Qn(f,f)(s,y)dy, t> 0, (4.1) 


in which t 0 is an arbitrary positive real number, and g 0 is a given function that 
stands for the initial data at time t a \ for practical reasons,we study GU rather 
than m that is the particular case when t 0 = 0 and g 0 = f 0 - Next, we prove, 
by a convenient choice of t Q that the global solution /„ of ED converges in a 
sense to specify, to a global solution / of 12.2011 . 

Proposition 4.1 Let r be an arbitrary strictly positive real number. 

Let g 0 £ X r and let n be an integer such that n > n 0 (r). 

Then, for every t 0 £ [0,+oo[, the integral equation has a unique solution 

f n £ C[t Q , +oo; X r ], Moreover, /„ satisfies the inequality : 



(4.2) 
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Proof: 

We give the proof in 2 steps. 

Step 1 :Local existence and estimation. 

Consider equation HD for t £ [t 0 ,t 0 + <5] where <5 > 0 is given. One verifies 
easily, using (3.8) that gives : Q n (f , /) + nf = n^/nR(n, Q)f and 


d \ £3 nt 
dt 


[e nt f{t 0 +1, y)\ = e r ‘ 


nf (t a + t,y), that (14.11) is equivalent to : 


f(to + t,y) = e nt g 0 (y)+ [ ny/ne n(t s) R(n, Q)f(t 0 + s, y) dy (4.3) 

Jo 

We solve ED by the fixed point theorem. Let us define the operator A over 
the complete metric space C[t 0 ,t 0 + 6;X r ] by : 

Af(t 0 + t,y) = e~ nt g Q (y) + f ny/ne~ n ^~ s) R(n,Q)f(t 0 + s,y) dy (4.4) 

Jo 

i) Let / £ C[t 0 ,t 0 + 5]X r \\ (14.411 gives, since || g 0 ||< r and using (3.9) that 
gives, || y/nR(n, Q)f(t a + s) ||=|| f(t 0 + s) || < r, and for every t £ [0, (5] : 

II Af(t a + t) || < e~ nt r + nre~ nt [ e ns ds = e~ nt r + re~ nt (e nt - 1) = r 

Jo 

Hence 11| Af \ ||< r, and this shows that A maps C[t a , t 0 + S ; X r \ into itself. 

ii) Let f,g £ C[t 0 ,t 0 + 6; X r ]; t£ [0,(5]; (4.4)gives: 

(Af - Ag)(t 0 +t,y) = [ ny/ne~ n ^~ s) [R(n, Q)f - R(n, Q)g] (t 0 +s, y) dy 

Jo 

which gives, using property (3.12) of R(n,Q) and e -n ( t-s ) < 1 : 


III Af - Ag |||< 2 nS ||| f-g ||| . 


It then appears that A is a contracting map in a metric space C[t 0 , t 0 + 5; X r \ 
where 2n6 < i.e. 6 £]0, j-]. Taking 6 = j-, we conclude that A has a 
unique fixed point f° £ C[t 0 ,t 0 + j^;X r ], that is the unique solution of (4.3) 
and hence, of (4.1). 

Now we have, by taking t = 0 in (4.3), f(t 0 ,y) = g 0 (y )■ The solution f° then 
satisfies: 


fn(to+t,y) = e nt f°(to,y) 



ns/ne n(t s) R(n,Q)f°(t 0 + s,y) dy 


which gives, multiplying by e nt , using once more (3.10), and for t £]0, 


e nt f° n (to + t) || <||/“|| +n 


e ns fn(to + S ) II ds 
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this gives, by Gronwall Lemma: e nt || /°(i 0 + t) ||< e nt || /°(f Q ) || and hence 


/nlll<ll/n(0 II 


(4.5) 


Step 2: Global existence; estimation and uniqueness. 

Let k G N. By replacing in the integral equation (4.1) t a by : 

t 0 + G + " • , to + , step 1 tells us that on each interval 

h = [to + t 0 + whose length is the initial value problem for the 

corresponding integral equation has a unique solution /* G C[/fc, X r ], provided 

that the initial data we denote /*; [t c + is a given element of f* then 

satisfies: 

j fn (to + it + t, y) = fh(to + ^,y) + ft°l¥ +t Qn (. fn > fn) 0 : 2 /) 
l/n(^ + |r)G^; fcGN; te[o,^] 

and (4.5)implies: 

\\\f k n \\\<\\f k n {to+^) II . 

Now, [0, +oo[= |J [t 0 + t 0 + £ + t] ; Define : 
keN 


fn (to +^ y) = fn 1 (to + y) if k > 1 and fn(to) = 9o- 

Then, the solution /^ _1 and that are defined on \t Q + , t a + -^] and 

[to + it'to + h tr) overlap at t, = t 0 +±. 

Define the function f n on [0, +oo[ by 


fn{t) 


fn(t) 


if 


t G 


to T . , t Q 
4 n 


k + 1 

4n 


Then, a straightforward calculation shows, using the above relations for 

k, k — 1, k — 2,.» ■ ■ ,1,0, that f n is a global solution of (4.1) on [0, +oo[ satisfying 

the estimation (4.2) and hence, f n G C[t 0 , +oo; X r ], 

Now suppose that f,g G C[t a , +oo; X r ] are 2 solutions of (4.1), then: 


(/ - 9 ) (to + t,y)= ( {Qn(f , /) - Qn(g, g)} (to + s,y)ds t> 0. 

J O 

This gives, using property (3.14) of Q n - 

II (f~g)(to + t) || < K ( || (/ - g) (t 0 + s) || ds t> 0 
Jo 

which implies, using Gronwall Lemma, / = g, and the uniqueness follows. This 
ends the proof of proposition 4.1.D 

We now state and prove the global existence theorem. 
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Theorem 4.1 Let f 0 € T 1 (M 3 ), f Q > 0, a.e; be arbitrary given. Then the 
Cauchy problem for the relativistic Boltzmann equation on the Bianchi type I 
space-time, with initial data f 0 , has a unique global solution f £ C[0, +oo; L 1 (K 3 )], 
/(f) > 0, a.e,'it £ [0, +oo[; / satisfies the estimation: 

Sup || f{t) || < || f 0 || (4.6) 

££[ 0,+oo[ 


Proof: 

Fix a real number r >|| f 0 ||. We give the proof in 2 steps. 

Step 1: Local existence and estimation. 

Let t 0 > 0 and g 0 £ X r be given. The proposition 4.1 gives for every integer 
n > n 0 (r), the existence of a solution /„ £ C[t 0 ,+ oo;X r ] of (4.1) that satisfies 
the estimation (4.2). It is important to notice that this solution depends on n. 
Let T > 0 be given; we have f n £ C[t 0 ,t 0 + T;X r \, Vn £ N, n > n 0 (r). Let 
us prove that the sequence f n converges in C[t 0 , t 0 + T] X r \ to a solution of the 
integral equation : 


f(t 0 + t,y) 


9o(y)+ f ~^Q (f, /) (s, y) dy, 


t £ [0, T] : 


(4.7) 


Consider 2 integers m,n> n Q {r) ; we deduce from (4.1), that, Vf £ [0, T]; 


fn(t 0 T t,y) ,fm{to T t, y) 



[Qn {fn,fn) 


Qm (fmi fm)] 2/) dy 


this gives, using the property (3.16) of the operators Q n , Qm and Gronwall 
lemma: 

III fn fm 11| — (- + -)KTe KT 
\n m J 

this proves that /„ is a Cauchy sequence in the complete metric space 
C[to, t 0 + T , X r ] . Then, there exists f £ C[t 0 ,t 0 + T, X r ] such that: 


f n converges to / in C[t 0 ,t 0 + T, X r ] 


(4.8) 


The definition of X r implies that f(t a + t) > 0 a.e Vt £ [0,T]; 

Let us show that / satisfies (4.7). The convergence (4.8) implies that, f n {t 0 + t) 
converges to f(t 0 +t), Vt £ [0, T]; then f n (t 0 ) converges to /(t D ); but f n (t 0 ) = g a , 
then /(t 0 ) = go', next we have, using property (3.14) of Q n and Q : 


nt 0 -\-t 


Qn {fn, fn) - Q (/, /) 

Vo 


(s, y) dy II < ^ + KT III /„ - f 


which shows that 

r t 0 +t 


/ t 0 ~\~t rt 0 -\-t 

Qn{fn,fn)(s)ds converges to j Q(f,f)(s)ds in L 1 (M 3 ). 
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Hence, since /„ satisfies (4.1), / satisfies (4.7), Vi £ [0,T] and a.e. in y £ R 3 . 
Finally, (4.2) and (4.8) imply: 


III /III <11 fito) II 


(4.9) 


Step 2 : Global existence, estimation and uniqueness 

Since t a £ [0, +oo[ is arbitrary and since Vn > n 0 (r), the solution f n of (4.1) 
is global on [ t a , +oo[, Step 1 tells us that, given T > 0, by taking in the integral 
equation (4.7), t 0 = 0, t a = T, t 0 = 2T, • • • ,t a = kT, • • • k £ N, then, on each 
interval Jk = [kT , ( k + 1)T], k £ N, whose length is T, the initial value problem 
for the corresponding integral equation has a solution f k £ C[Jk, A' r ], provided 
that the initial data we denote f k (kT), is a given element of X r . 

We then proceed exactly as in Step 2 of the proof of Proposition 4.1 by writing 
for T > 0 given, that [0, +oo[= (J [kT, (k + 1)T], to overlap the local solutions 


f k £ C[kT, (k+l)T;X r \, by setting f k (kT ) = f k ~ 1 (kT) if k > 1 and f°( 0) = f a , 
and we obtain a global solution / £ C[0, +oo; X r ] of the integral equation (2.20) 
defined by f(t) = f k (t), if t £ [kT, (k + 1)T]. Next, (4.9) that gives 


III f k III < II f k {kT) || 


shows that / satisfies the estimation (4.6). 

Now if /, g £ C[0, +oo; L 1 (R 3 )] are 2 solutions of (2.20) for some initial data f a , 
(4.6) and property (3.5) of Q give: 



which gives, by Gronwall Lemma, f = g, and the uniqueness follows. 

Now as we indicated, solving the initial value problem for the relativistic Boltz¬ 
mann equation is equivalent to solving the integral equation (2.20). This ends 
the proof of theorem 4.ID 

Remark 4.1 1) The expression C(p,q,uj) = u>.(p — q ) used by the author 

in m is not correct in the full relativistic case, and is to be replaced by 
(2.15) in which the dot product that corresponds to the Bianchi type I case 
is given by (2.3) . 

2) In ITU( the author defined the operator R(n , Q ) whose domain is X r , using 
the equation: nu — A- Q(u,u ) = v, instead of (3.6), and he associated to 
it, the approximating operator P n {u,u) = nR{n,Q)nu —nu, instead of Q n 
defined by (3.8). But it is clear that, even if u £ X r , nu is no longer in 
X r for n sufficiently large. So the domain of P n is unspecified. 

3) One could use (3.5) that means that Q is locally lipschitzian in f, to solve 
the last equation (2.17) using classical results on the differential equation. 
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But this method does not give the positivity property of f, imposed by 
the physical nature of the considered problem, and that we obtain by our 
method. 

4) In theorem f.l, there is no restriction on the size of f 0 which can be taken 
arbitrary large in L 1 — norm 
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